Test bodies and naked singularities: is the self-force the cosmic censor? 
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Jacobson and Sotiriou showed that rotating black holes could be spun-up past the extremal limit 
by the capture of non-spinning test bodies, if one neglects radiative and self-force effects. This would 
represent a violation of the Cosmic Censorship Conjecture in four-dimensional, asymptotically flat 
spacetimes. We show that for some of the trajectories giving rise to naked singularities, radiative 
effects can be neglected. However, for these orbits the conservative self-force is important, and 
seems to have the right sign to prevent the formation of naked singularities. 



PACS numbers: 04.70.Bw, 04.20.Dw 

The most general stationary vacuum black-hole (BH) 
solution of Einstein's equations in a four-dimensional, 
asymptotically flat spacetime is the Kerr geometry 
characterized only by its mass M and angular momentum 
J. Solutions spinning below the Kerr bound cJ/GM 2 < 
1 possess an event horizon and are known as Kerr BHs. 
Solutions spinning faster than the Kerr bound describe 
a "naked singularity", where classical General Relativ- 
ity breaks down and (unknown) quantum gravity effects 
take over. It was hypothesized by Penrose that classical 
General Relativity encodes in its equations a mechanism 
to save it from the breakdown of predictability. This is 
known as the Cosmic Censorship Conjecture (CCC) Q, 
which asserts that every singularity is cloaked behind an 
event horizon, from which no information can escape. 

There is no proof of the CCC. Indeed there are a 
few known counter-examples, but these require either ex- 
treme fine-tuning in the initial conditions or unphysical 
equations of state Q , or are staged in higher-dimensional 
spacetimes 0. Moreover, all existing evidence indicates 
that Kerr BHs are perturbatively stable [|J, while Kerr 
solutions with cJ/GM 2 > 1 are unstable Q. Thus, 
naked singularities cannot form from BH instabilities. 

Because naked singularities appear when cJ/GM 2 > 
1, it is conceivably possible to form them by throwing 
matter with sufficiently large angular momentum into a 
BH. With numerical-relativity simulations, Ref. @ found 
no evidence of formation of naked singularities in a high- 
energy collision between two comparable-mass BHs: ei- 
ther the full nonlinear equations make the system radiate 
enough angular momentum to form a single BH, or the 
BHs simply scatter. The case of a test-particle plunging 
into an extremal Kerr BH was studied by Wald 0, who 
showed that naked singularities can never be produced, 
because particles carrying dangerously large angular mo- 
mentum are just not captured. 

Recently, Jacobson and Sotiriou (JS) @| (building on 
Refs. 0]) have shown that if one considers an almost ex- 
tremal BH, non-spinning particles carrying enough angu- 
lar momentum to create naked singularities are allowed 



to be captured. 1 As acknowledged by JS, however, their 
analysis neglects the conservative and dissipative self- 
force (SF), and both effects may be important In 
this letter we will show that the dissipative SF (equiv- 
alent to radiation reaction, i.e. the energy and angular 
momentum losses through gravitational waves) can pre- 
vent the formation of naked singularities only for some 
of JS's orbits. However, we will show that for all these 
orbits the conservative SF is comparable to the terms 
giving rise to naked singularities, and should therefore 
be taken into account. Hereafter we set G = c = M = 1. 

JS considered a BH with spin a = J/M 2 = 1 — 2e 2 , 
with e -C 1, and a non-spinning test-particle with energy 
E, angular momentum L and mass m. Neglecting the 
dissipative and conservative SF, the particle moves on a 
geodesic, and JS identified a class of equatorial geodesic 
orbits such that (i) the particle falls into the BH, which 
implies an upper limit on the angular momentum, L < 
L mm , and (ii) the BH is spun up past the extremal limit 
and destroyed, which implies a lower limit on the angular 
momentum, L > L mi „. Therefore 
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2E + E < L < L T 



(2 + 4e)E. (1) 



Imposing L max > L m i n then yields 

E min = (2 - V2)e <E< E max = (2 + y/2)e . (2) 

Finally, JS checked that these intervals contain both 
bound orbits (i.e. orbits that start with zero radial veloc- 
ity at finite radius) and unbound orbits (i.e. orbits that 
start from infinity) . Parameterizing the above interval as 



E = E, 
L = L„ 
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(3) 



= imin + 8ye 2 (l - x)x (4) 



JS consider also the case of spinning particles, but in this letter 
we will focus on the non-spinning case. 
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with 0<x<l,0<y<l, the final spin is 



n JS _ 
a f ~ 



L 



(1 + Ef 



1 + 8e 2 (l -x)xy + 0(e 3 ) > 1, (5) 



and the spin-up is due to the terms quadratic in e. 

Let us first investigate how radiation reaction affects 
JS's analysis. Taking radiation losses -E ra d an d -£<rad into 
account, Eq. becomes 

a f = 1 + 8e 2 (l - x)xy + 2E Iad - L rad + 0(e 3 ) . (6) 

Let us focus on unbound geodesies, 2 and following JS 
assume E/m^> 1 and L/m^$> 1 (null orbits). These or- 
bits are characterized by the impact parameter b = L/E 
alone. From Eqs. ([1]) and ©, JS's orbits have L = bE, 
with b = 2+4e[l-2x(x-l)(j/-l)/(2+V5(2x-l))]. Vary- 
ing x and y between and 1, one obtains b = 2 + 5e, with 
2y/2 < 5 < 4. However, because b ph = 2 + 2y/3e+(D(e 2 ) is 
the impact parameter of the circular photon orbit ( "light- 
ring"), only orbits with 2v2 < 5 < 2y/3 are unbound. 
When S as 2^3, these orbits are expected to circle many 
times around the light-ring, so radiation reaction could 
prevent the formation of naked singularities or at least 
invalidate JS's analysis. In fact, for 5 arbitrarily close to 
2-\/3, the particles would orbit around the light-ring an 
arbitrarily large number of times, and gravitational- wave 
emission must be important (Tlj . We will show, however, 
that this is not true for all of JS's orbits. 

Considering the geodesic equations for null equatorial 
orbits with impact parameter b = 6 pri (l — fc), with k <C 
e -C 1, one finds that the radial potential - defined as 
V r (r) = (dr/dX) 2 with A an affine parameter - has a 
minimum at r = r imn = r p h + O(k), near which 




—fee + 3(r - r min ) 2 
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Integrating from r min — Ar 2 to r m - m + Ar±, with Ar\ >2 3> 
ke, the number of cycles near the minimum is 



N, 



cycles 



r min +Ari 



-Ar 2 



-- = [^+51o g (*e)] - + - 



(8) 

A and B being constants depending on the integration 
interval. Fixing e, and thus the BH spin, we can see that 
-^cycles depends on logfc, and diverges when k — > [Tlj . 
Because the fluxes are proportional to 7V cyc i 0S , we have 



£ rad = AE(e) x N, 



cycles i 



L 



rad 



AL(e) x iV cyclcs , (9) 



2 As we mentioned, JS also considered bound orbits, falling into 
the BH from a Boycr-Lindquist radius r = r^ OI +0(e) (r hor being 
the horizon's radius). However, these orbits pose a problem, as 
we will show later, because the distance to the horizon is compa- 
rable to the particle's minimum attainable size max(E, m) > e, 
so finite-size effects should be taken into account. 



where AE and AL are the fluxes in a single orbit. From a 
frequency-domain analysis \12\ . AE/AL must equal the 
light-ring frequency, il p h m 1/2 — (V3/2)e, hence 



AE{e) = Ex{e){l + e 2 e), 
AL(e) = 2E 1 (e)[l + (V3 ■ 



e2)e] 



(10) 
(11) 



Here -Ei(e) is the energy flux for a single orbit at leading 
order in e, and e 2 is an undetermined coefficient. Semi- 
quantitative arguments by Chrzanowski [13| and more 
rigorous analytical calculations by Chrzanowski and Mis- 
ner fijl show that, 



Ei - (r-r H )E 2 - eE 2 - e 3 



(12) 



(later we will discuss an additional proof of this scaling). 
At leading order in e, this results in 



£ rad = AE(e) x N cyclcs - log (ke)e 2 



(13) 



This scaling still depends on k, but the dependence is log- 
arithmic, so unless k is really small -E ra d ~ log (e)e 2 . Al- 
though terms of order e 2 log e seem to dominate Eq. © , 
because of Eqs. ©, JIDJ and (JTTJ> one has i ra d 2i? rad — 
2\/3ei?i(e)-/V cyc i cs ~ e 3 loge. Therefore JS's analysis is 
valid for these trajectories. However, if k < exp(— 1/e), 
-E"r a d ~ e and ^ ra d — 2i? ra( j ~ e 2 , and JS's analysis is not 
valid because radiative effects cannot be neglected. 

To test the above picture we used a time-domain 
code 11511 solving the inhomogencous Teukolsky equa- 
tion [l2| that describes the gravitational perturbations 
of Kerr BHs in the context of extreme mass-ratio bina- 
ries (EMRBs). This code has been successfully used in 
many scenarios, inclu ding an extensive study of recoil ve- 
locities from EMRBs [Ty| . Because, for almost extremal 
BHs and in Boycr-Lindquist coordinates, the particle's 
orbit, the light-ring and the horizon are extremely close, 
we modeled the test-particle to have a fixed width in the 
"tortoise" coordinate r* as opposed to r [13] , and checked 
that our results are independent of the particle's width 
when that is sufficiently small. (More details on these 
tests will be presented in a follow-up paper.) 

We consider BHs with a = 0.99, 0.992, 0.994, 0.996, 
0.998 and 0.999 and geodesies having E = (E max + 
£min)/2 = 2e, L = b ph E(l - k) with k = 10~ 5 , and m = 
0.001 <C E. Using these geodesies and integrating their 
cycles from r = 1.05r ph to r = (r ph + r hoj: )/2 (r ph being 
the light-ring radius), we get A 0.3294, B m -0.01941 
for the coefficients in Eq. ([SJ. Assuming E\ = e-ye n , we 
fit the energy and angular-momentum fluxes at infinity 
with Eqs. (|9^)- pT|) . obtaining n xs 2.91. Because this is 
very close to the theoretical value n — 3, we assume n = 3 
and fit the data with only two free parameters, e\ and e 2 , 
obtaining ei = 136.97 and e 2 = —4.423. With these val- 
ues, Eqs. (|9^- (fl"Tj) reproduce the numerical data to within 
1 — 3% for a < 0.999, which is comparable to the data 
accuracy. For a = 0.999, however, the fluxes predicted 
by Eqs. ©-fll]) arc about 12% larger than the numeri- 
cal ones. To investigate this issue, we ran an additional 
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TABLE I: Initial and final BH spin after absorbing a par- 
ticle with energy E — y2(l ~~ a ) an< A angular momentum 
L — b p hE(l — 10 _J ), neglecting conservative SF effects, but 
not radiation reaction. We also show the final spin without 
radiation reaction (aj s ) predicted by JS. 
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0.99 


0.992 


0.994 


0.996 


0.998 


0.999 


0.9998 


«/ 


0.882 


0.928 


0.961 


0.984 


0.997 


0.9996 


1.00006 




1.0043 


1.0035 


1.0026 


I.00I8 


1.0009 


1.00045 


1.00009 



simulation for a = 0.9998, which seems to confirm that 
Eqs. (|9"fl- (fTTj) ovcrprcdict the fluxes for very high spins. 
At this stage it is not clear whether this is a numerical 
problem (simulations are very challenging for a w 1) or 
whether this is due to the simplified analytical deriva- 
tion of Eqs. (|fI|)- ([TT|) . We will investigate this issue in the 
follow-up paper, but because the numerical fluxes are 
smaller than expected, it only reinforces our conclusion 
that there are orbits giving rise to naked singularities 
even when radiation reaction is taken into account. 

Since L rad - 2£ , rad = 2V3eE 1 (e)N cyclcs ~ e 3 log(fce) > 
0, Eq. predicts that radiation reaction will decrease 
the final spin a/. Using the above values for A, B 
and ei, and x = 0.5 and y ps 2-^/3 — 3 + 4e/3 corre- 
sponding to our geodesies, Eq. (|6|) predicts a/ < 1 for 
e ^ e crit ~ 0.003. However, for sufficiently large spins, 
the term L rad — 2£ , ra( j ~ e 3 log(e) is subdominant and 
a/ > 1. Numerical results confirm this expectation: in 
Table HI we show the BH spin a/ after absorbing the par- 
ticle, taking into account radiation reaction. As can be 
seen, a/ > 1 already for a = 0.9998, corresponding to 
e = 0.01 > e C rit- This is because, as already mentioned, 
Eqs. ([9])- ([IT]) overpredict the fluxes for a > 0.999. 

Moreover, even the result that af < 1 for e > e cr jt 
is questionable. Indeed, the fluxes down the horizon 
might destroy the BH before the particle is captured, 
while our code only calculates the fluxes at infinity. This 
is sufficient for our purposes because we used the code 
only to test the scaling (|9))- (|TT|) . which is expected to 
hold both for the fluxes at infinity and down the hori- 
zon, since its derivation is generic. Once validated, that 
scaling implies that for sufficiently large spins both fluxes 
are smaller than the terms giving rise to naked singular- 
ities [i.e. the quadratic terms in Eq. ||BJ)]. For e > e cr i t , 
instead, the fluxes at infinity decrease the final spin to 
dfin < 1. However, in such a situation also the fluxes 
down the horizon, I/ ra d,in and -Erad.in, are expected to 
be important (because the fluxes are produced when the 
particle sits at the light-ring, which roughly corresponds 
to the maximum of the effective potential for gravita- 
tional waves), and could destroy the horizon before the 
particle is captured. In fact, the spin change is Aa = 

£rad,in — 2£/ ra d,i n = £'rad,in(l/^ph — 2) PS 2-\/3e-Erad,i ri , be- 
cause -Eradun/irad.in = O ph PS 1/2 - (\/3/2)e. Since fi ph 

is larger than the horizon's frequency flh m ~ 1/2 — e, 
radiative emission is non-supcrradiant and £rad,in > 0, 
hence Aa > 0. Thus, the ingoing fluxes increase af. 



So far we have shown that radiation reaction cannot 
prevent the formation of naked singularities, unless the 
impact parameter is extremely close to the light-ring's 
impact parameter & pn . We will now show, however, that 
for all of JS's orbits the conservative SF is as important 
as the terms giving rise to naked singularities. 

Let us consider a BH 3 with gravitational radius R g = 
2Gm/c 2 in a curved background with curvature radius 
£ ^> R g . The rigorous way of studying the motion of 
this BH is to set up a proper initial value formulation, but 
a reasonable alternative for practical purposes is to use 
a matched asymptotic expansion [l8| . Near the BH (i.e. 
for r < Ti, ri being a radius <C £), the metric is (/internal = 
5bh + H\(r/C) + H 2 (r/ C) 2 + where <7bh is the metric 
of an isolated BH and Hi(r/£), H2(r/£) 2 are corrections 
due to the "external" background. Far from the BH (i.e. 
for r > r e , r e being a radius 3> Rg), the metric is in- 
stead ^external = ^background + ^1 (Rg I £) + ^2 (Rg / £) + ■ ■ • ! 

i.e. the background metric plus perturbations due to the 
BH's presence. Because R g <C C, there exists a region 
r e < r < ri where both pictures are valid and the two 
metrics can be matched. Doing so, one finds that the BH 
equations of motion are [l8| 

«"V/ = /cons + /dk + 0(R g /C) 2 , (14) 

where V is the connection of the background spacetime. 
The terms /^ ons and f^ ss are G(R g /C), and are known as 
the conservative and dissipativc SF. Remarkably, it turns 
out that Eq. (fT4"|) is the geodesic equation of a particle in 
a "perturbed" metric g = g + h R , where h R is a smooth 
tensor field of order 0(R g /C): 

u^V^u" = (15) 

(the connection V and the 4-velocity u M being defined 
with respect to the "perturbed" metric g = g + h R ). 

The dissipative SF amounts to the energy and angular- 
momentum fluxes considered earlier. Taking for instance 
the energy loss, Eq. (fl"4|) and E = —pt give dE/dr = 
— m/ t dlss = 0(R g /L) 2 . Assuming now that the back- 
ground spacetime is a BH with mass M ~ £ ^> R g , and 
specializing to orbits near the horizon, one has dt/dr ~ 
r H /(r - r H ), which gives dE/dt ~ (r - r H )0(R g /C) 2 . 
Comparison of this scaling with our numerically validated 
scaling (fT2"|) shows that for a BH with E 3> m the size 
entering the matched asymptotic expansion above (the 
"physical" size) is R g = 2GE/c 2 and not R g = 2Gm/c 2 . 
This is no surprise, as the physical size associated with an 
ultrarelativistic BH is dictated by its energy and not by 
its mass, because in General Relativity energy gravitates. 
Remarkably, however, we were able to test this fact with 
the numerical results presented earlier. Further evidence 



This discussion is completely general because the motion of a 
BH is the same as that of a particle with mass m, at leading and 
next- to-leading order in R g /C [T3l . 
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comes from boosting the Schwarzschild line-element to 
the speed of light, keeping the total energy fixed. One 
gets the Aichelburg-Sexl metric, which depends on the to- 
tal energy E and not on the rest-mass [19| : this boosted 
BH absorbs particles within a distance ~ E from it. 

Because a BH's size is determined by max(E,m) > e, 
the conservative SF affects JS's analysis. This is easier 
to see from Eq. (fT5j) (although the same result can be ob- 
tained from Eq. (fl4|) : see Ref. [2(| ): because the metric 
"perturbation" h R is 0(R g /C) = 0(e), the effective po- 
tential for the radial motion differs from the "geodetic" 
one by 0(R g /C) = 0(e) [HS2. Therefore, b ph changes 
by Sb p h = G(R g /C) = 0(e). Because JS's orbits have 
frph — b — 0(R g /C) = 0(e), the conservative SF may 
prevent them from plunging into the horizon. This effect 
is intuitive: if the particle's size is ~ e, finite-size effects 
are important for impact parameters b = b p h + 0(e). 

A calculation of Sb p h is not doable with present tech- 
nology [22|, but we can estimate its sign. Because of 
frame-dragging, for a = 1 — 2e 2 one has o p h = l/fi p h + 
0(e) 2 . While the SF effect on f2 p h has not been calcu- 
lated yet, Ref. [2(| calculated the Innermost Stable Circu- 
lar Orbit (ISCO)-frequency shift for a = 0, and showed 
that the conservative SF increases fiisco- It therefore 
seems plausible that f2 p h should follow the same behav- 
ior. While approximate methods for calculating the con- 
servative SF in Kerr spacetimes exist 0, [24| , they have 
problems for large spins, and the definitive answer to 



whether 6 p h increases or decreases for awl will only 
be available when a rigorous SF calculation [22[ is per- 
formed. However, assuming that the a = behaviour of 
Ref. [2Q] holds also for awl, one obtains that Q p h in- 
creases due to the SF, and therefore o p h should decrease, 
possibly preventing the capture of the particles with dan- 
gerously large L and the formation of naked singularities. 

In conclusion, we have shown that radiation reaction 
effects can prevent the formation of naked singularities 
only for some of the orbits for non-spinning particles 
around almost extremal Kerr BHs identified by JS. How- 
ever, for all orbits capable of producing a naked singu- 
larities, the conservative SF is non-negligible and seems 
to have the right sign to prevent the particles from being 
captured, thus saving the Cosmic Censorship Conjecture. 
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